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We  have  shown  recently  that  the  steady- state  behavior  of  simple  reactions 
at  microdisk  electrodes  can  be  predicted  '  by  the  application  of  Neumann's 
integral  theorem(3>.  The  concentration  distribution  for  a  reactant  (or 
product)  in  such  a  reaction  is  governed  by 


Di!£  +  Dac  +  Da^c_0 

dr2  r  dr  dz2 


and  we  obtain 


C(r ,z) 


.  f.d.r*. 

*  ft  *  A 


exp(-az)  Q(p) 


Jo(ar)Jo(ap)pdp 


Equation  [2]  makes  a  formal  linkage  between  the  distribution  of  sources,  Q(p) 
(or  sinks)  at  the  radial  position  p  and  the  concentration  throughout  space. 

We  have  shown  that  the  use  of  the  simple  form  for  the  distribution  of  sinks 


-Q(tf)  - 


a  cos(0) 


—  { 
v 


c  +  c  cos (8)  +  c  cos(20)  +  •••  +  c  cos(n0) 

0  1  2  n 


where 


p  -  a  sin(0) 


2 


leads  to  the  recovery  of  known  results(<l'5>(derived  by  the  application  of  the 
discontinuous  integrals  of  Bessel  functions)  for  constant  surface 
concentrations  provided  we  assume 

c  -  c  -  c  -  •  •  •  -  c  -  0  [5] 

12  3  n 


and 


and  of  constant,  uniform,  surface  flux  provided  ve  assume 


c  —  c  —  c  —  •••—  c  —  0 

0  2  3  n 


[6] 


[7] 


and 


c  -  Qa 

l 


[8] 


It  has  also  been  shown  that  the  reexpression  of  [3]  in  the  form 


-Q(*>  - 


a  cos(0) 


d  +  d  co s(8)  +  •••  +  d  cos  (0) 
0  1  n 


[9] 


is  particularly  useful  in  discussing  the  application  of  more  general  boundary 


conditions  such  as 


D(— )  -  kcS<  0  <  r  <  a,  z-0 
KdzJ 


[10] 


The  use  of  [9]  in  [2]  leads  to  the  simple  result 


C(r) 


i/2  V"  r~  r,  . 

-  -  - -  )  d  -4-^-  pft.-i;  1;  -  \ 

2D  /  J  2  H2  2  a2/ 

1 I  2 

J-0  V 


[11] 


(at  z-0)  where  T  is  the  gamma  function  and  denotes  the  appropriate 
hypergeometric  function.  Application  of  [11]  at  (n+1)  radial  positions  in  any 
particular  boundary  condition  (such  as  [10])  then  allows  the  evaluation  of  the 
unknown  coefficients  in  [9]  thereby  giving  a  complete  solution  for  any 
particular  problem.  The  results  obtained  by  the  application  of  Neumann's 
integral  theorem  has  been  related  also  to  results  derived  by  an  extension<6)of 
the  application  of  the  discontinuous  integrals  of  Bessel  functions**’55 .  We 
have  shown  also  that  the  use  of  Neumann's  integral  allows  the  evaluation  of 
the  effects  of  the  distribution  of  potential  in  the  solution  (the  tertiary 

current  distribution)  and  of  the  non-linearities  in  the  concentration  terms  o : 

(1  2) 

the  electrode  reactions 

In  this  paper  we  extend  the  approach  to  the  discussion  of  the 
chronoamperoraetric  response  of  microdisk  electrodes,  comment  on  the  effects  o: 
the  distribution  of  potential  in  the  solution,  and  outline  the  analysis  of 
linear  sweep  voltammetric  experiments.  In  view  of  the  formal  similarity  of 
the  derived  equations  to  those  obtained  for  the  steady  state  behavior  of 
reactions  following  the  c.e.  and  e.c.  mechanisms,  we  also  derive  the  relevant 
results  for  these  reaction  pathways. 


We  consider  first  of  all  the  behavior  of  a  simple  diffusion  controlled 


reaction  where  mass  transfer  of  the  reactant  is  controlled  by 


ac 

D  — 

at 

ar2 

r 

ar 

az2 

[12] 


subject  to  the  initial  condition 

C-C®,  0  <  r  <  ®,  z  >  0,  t-0 

and  the  boundary  conditions 

C-0,  0<r<a,  z  -  0,  t>0 

3C 

0  — ■—  —  0  z  .  o  t  >  0 

dz 

Laplace  transformation  of  [12]  subject  to  [13]  gives 


[13] 


[14] 


[15] 


a2c  |  l  ac  |  azc  s  -  |  c" 

drz  r  dv  dz2  D  D 


0 


[16] 


where  s  is  the  variable  of  the  Laplace  transformation.  The  solution  of  [16] 
can  therefore  be  expressed  in  terms  of  Neumann's  integral  of  two  parameters 


i 


5 


(p,s)  using  Q(p,s)  as 


1  C(r,z,s) 

By  analogy  with  the  expansion  [9]  used  to  describe  the  distribution  of  sinks 
I  in  the  steady  state,  ve  write 

r'Q <*)  - - - 

a  cos (6  ) 


d  (s)  +  d  (s)cos(0)  +  ••• 

0  X 


+  d  (s)cos  ( 8 ) 
n 


[18] 


oo 

-  . 


a  da 


f  [  2  si 

Tla  +  5|  z[ 

[1  1/2 

”!+l] 


Q(p,s) 


JQ(ar)Jo(ap)pdp 


[17] 


► 


where  all  the  coefficients  d  are  now  functions  of  s.  Substitution  of  [18 

J  1 

and  [4]  in  [17]  gives  at  z  -  0 


00 

z,  /  \  C  a 

C(r , s)  -  - 


Jq (ar)ada 
1/2 


r*72 


d  j  r  2  si 

0  L“ +  sj 


n 

I*. 


j-0 


(s)J  (aa  sin(0))sin(0)cos  ($)d8 
o 


6 


C® 


s 


n 


J-0 


J  (ar) 
0 


(aa)  d(aa) 


( ( j -l > / 2  > 


[19] 


for  a  distribution  of  sinks  Q(0)  over  the  surface  where 


e^  -  2<<J'1>/2)r((J+l)/2) 


[20] 


The  coefficients  d^(s)  can  then  be  obtained  by  substituting  [18]  and  [19]  into 
the  Laplace  transform  of  the  boundary  condition  appropriate  to  the  particular 
form  of  the  experiment.  For  example,  for  the  simple  condition 


C-0,  0  <  r  <  a,  z-0  [21] 


(cf  [14])  we  solve  the  set  of  equations 


n 


J  (ar  )  J, 
o  k  f 


J  +  i'l 


(aa) 


(aa)  d(aa) 


(  ( j - 1) / 2  ) 


1 


s 


[22] 
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at  (n+1)  radial  positions.  We  therefore  obtain  the  unknown  coefficients  as 
functions  of  the  parameter  s.  Substitution  of  these  coefficients  into  [18] 
and  integration  over  the  radius  of  the  disk  then  gives  the  Laplace  transform 
of  the  total  rate  of  reaction,  R, 


n 


J-0 


(s)cosJ 


(Osin(tf)d* 


T  Vs) 

2wa  1  - 

zL  u+n 


[23] 


Numerical  inversion  of  [23]  gives  the  total  rate  of  reaction  in  the  t-domain. 

An  alternative  procedure  is  to  invert  [22]  to  the  relevant  convolution 
integrals 


n 

I 


j*0 


A 


-.1/2  1/2 

D  n  a 


J  (or  )  J 

0  k 


o  J 


d  (t-r)  exp (-Da  r) 

„  J  1/2 
0  T 


J+l 

2 


(aa)  d(aa)dr 


faa) 


(( j -  1  )  12) 


d  (r)  exp (-Da  (t-r)) 


3  (t-r) 


J  (ar  )  J ,  \ (aa)  d(aa)dr 

o  k  I  J  +  l] 


( ( j  - 1 ) 12  ) 


In  the  numerical  integration  «ver  £ ay  1  intervals  of  Ar  we  then  have  to  solve 
the  set  of  Equations  [24]  for  the  values  of  d^  in  the  1th  interval.  This 
forward  integration  is  simplified  iry  the  fact  that  we  know  that  the  system 
follows  the  Cottrell  equation  at  short  times<7>for  which  at  sufficiently  large 


j  /  \ 

d  (s)  -  - 

1  1/2 
S 


d  —  d  —  d  —  •  • 

0  2  3 


-  d  -  0 


Other  boundary  conditions :  the  tertiary  current  distribution;  linear  sweep 
voltammetry . 

We  have  pointed  out  elsewhere  that  the  application  of  the  boundary 
condition  [14]  (or,  more  generally,  of  a  constant  surface  concentration 
boundary  condition)  is  unrealistic  as  the  flux  would  have  to  become  infinite 
at  the  edges  of  the  disk'1, 2,*'8) ; see  also<6) .  However,  the  flux  cannot 
become  infinite  for  two  reasons:  firstly  because  the  rate  constant  governir. 


9 


the  electrochemical  reaction  must  itself  be  finite  and,  secondly,  because  an 
|  infinite  rate  of  reaction  would  alsj  require  an  infinite  overpotential  at  he 

edges  of  the  disk.  The  current  distribution  must  therefore  be  more  uniform 
than  that  predicted  using  the  boundary  condition  [14]  under  most  experimental 
I  conditions  (the  tertiary  current  distribution). 

As  the  application  of  Neumann's  integral  theorem  provides  a  link  between 
the  assumed  form  of  the  flux  distribution,  Equation  [18],  and  the 
|  concentration  distribution,  Equation  [19],  it  is  straightforward  to  assess  the 

effects  of  the  boundary  conditions  more  complex  than  [14].  Ve  restict 
attention  here  to  irreversible  reactions  which  are  described  by  [10], 

•  Substitution  of  [19]  into  the  Laplace  transform  of  [10]  gives 


10 


s 


n 


J-o 


Jfl(ar)  J  /  +  x (aa)  d(aa) 

_ LID _ 

r  t  1 12 

I  2^  S  .  . ( ( J-l) 12) 

+  dJ  (aa) 


[27] 


The  dependence  o£  the  coefficients  d^  on  s  can  therefore  be  obtained  as  for 
^  the  case  of  zero  surface  concentration  governed  by  Equation  [22].  The 

distribution  of  the  flux  and  concentration  and  the  total  rate  of  reaction  can 
then  be  derived  for  the  Laplace  plane  using  Equations  [18],  [19],  and  [23]. 

Numerical  inversion  gives  the  dependencies  in  the  t-plane.  It  therefore 

I 

becomes  possible  to  assess  the  effects  of  changes  in  k  on  the  behavior  of  the 
system  and  the  relaxation  of  the  rate  of  reaction  under  potentiostatic 
^  conditions  can  be  evaluated. 

The  application  of  Neumann's  integral  theorem  also  allows  the  evaluation 
of  the  potential  distribution  in  the  solution,  as  in  the  case  of  the 
^  discussion  of  the  steady  state  However,  in  contrast  to  that  case,  the 

effects  of  this  potential  distribution  on  the  chronoamperometric  transient  can 
only  be  evaluated  in  limiting  cases.  The  potential  distribution  in  the 

presence  of  excess  support  electrolyte  is  governed  in  the  Laplace  space  by 

I 


d2j>  t  k  d]>  dZ4> 

K  -  +  -  —  +  K.  - -  - 

3r2  r  3r  dz2 


4>  -  0,  r 


«,  z 


GO 


[28] 


I 

while  the  distribution  of  sources  (for  a  cathodic  reaction)  over  the  disk  is 
given  by 

1 


» 

1 

I 

i 
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K 


di 


dz 


-zFQ(tf) 


-zF 

a  cos (8  ) 


d  (s)  +  d  (s)cos (0)  +  •••  +  d  (s)cos”(0) 

0  1  n 


[29] 


We  therefore  obtain  the  potential  at  the  surface  of  the  disk,  z  -  0,  as 


*<r,s) 


,s)  JQ(ar)  JQ(aa  sin(0)  ain(fi)  cos  8  &8da 


-zFa 


J  (arjda 
0 


II 


j-o 


r*/2 


d^(s)  JQ(aa  sin(0)  sin(0)  cos d 8 


n 


J-0 


where  e^  and  1^1  are  given  as  in  the  steady  state  case(1,by 


[31] 


-  2<<J'1,/2)r((j+l)/2) 


and 

I 


f 

J 


1/2 


2<<4*n/2>r((j+2)/2)a 


4? 


[32] 


Modification  of  [10]  to  take  Into  account  the  distribution  of  potential  in  the 
solution  gives 


.  k  exp[_^_(r,t)Flcs 
1  dz  >  RT  » 


(  <*</>( r,t)F  | 

V  RT  I 


-  k  exp|-“-^^  1(C~-  C(r,t)) 
RT 


[33] 


where  k  is  the  value  of  the  rate  constant  of  the  electrode  reaction  in  the 
absence  of  any  ohmic  potential  drops  in  the  solution  and  C(r,t)  is  the  inverse 
of  the  second  term  in  quation  [19],  A  general  Laplace  transform  of  [33]  can 
evidently  only  be  obtained  for  sufficiently  small  perturbations  of  the 
potential  and  concentration  for  which 


i 
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fac(r,t)'| 

a  k 

L  ,  o^(r,t)F  ] 

(C®-  C(r,t)) 

'  dz  ' 

RT 

a  k 

f(C*-  C(r , t) )]  + 

ka*(r,t)FC® 

l  ) 

RT 

[34] 


giving 


D 


1 

a  coa(9) 


d  <s)  +  d  (s)cos(#>  +  •••  +  d  (s)cos“(0) 

0  1  n 


fsa2] 

e 

i  D  ■ 

i 

>"»  JQ (or)  Jfx ^(aa)  d(aa) 


M 


J-0 


[•*♦  §] 


1/2 


(aa) 


(  U -i > / 2  ) 


2  00 

zaF  kC  a 
*RT 


n 


[35] 


j-0 


Special  solutions  using  [33]  include  those  cases  where  the  d^can  be 


approximated  by  linear  functions  in  time 


I 


dj  -  6j  +hJ<t) 


A  further  interesting  case  related  to  [36]  is  that  of  linear  sweep 
voltammetry .  We  again  restrict  attention  to  the  case  of  an  irreversible 
reaction  and,  neglecting  the  effects  of  the  tertiary  current  distribution, 
write  [10]  aa 


Df3C(£lt)j  -  k  exp[^](C*-  C(r,t)) 

V  t  V  RT  ) 


-  k  exp  (/it) (C  -  C(r,t)) 


where  k  is  the  rate  constant  at  t  -  0.  Laplace  transformation  of  [37]  gives 


D(-) - — 

Kdz'  a  cos (0) 


d  (s)  +  d  (s)cos(0)  +  •••  +  d  (s)cos  ( 6 ) 
0  1  n 


(s-/i) 


n 

— y~a 
Da  L J 


JQ (ar)  J  /  +  x (aa)  d(aa) 


(s-p)  e 


( < j  -  1  )  /2) 


[38] 


an  equation  system  which  has  to  be  solved  in  the  t-domain  where 


a  cos(ff) 


d  (t)  +  d  (t)cos(0)  +  •••  +  d  (t)cos  (0) 

0  1  n 


kC*exp  (/it)  -  — 


igw  ye  [[d 
>'V/2a  /  J  J  J 

‘ — -  •'o’'  o 

j-  0 


Jfl(or)J (aa)d(aa)dr 


(t-r)  exp( -Do  r ) 
_i/2 


2 


(aa) 


<<J-l>/2> 


-  kC 


'*XI,<''t)  •  ^TTFifr  y*i  f  k(,) 

»  *  *Z_  J„J0 

j-  o 


Jo(ar)J,,  ^ (aa)d(aa)dr 


O) 


f— 1 
2 


(aa) 


( ( j-l)/2) 


[39] 


c.e  and  e.c.  catalytic  reactions  in  the  steady  state 

We  first  consider  the  behavior  of  the  c.e.  reaction 


[i] 


ze 


[ii] 


in  the  steady  state  with  the  reactant  A  present  in  excess  concentration.  Then 
if  C  denotes  the  concentration  of  B,  the  diffusion  of  this  species  is  governed 


16 


by 


Di!c  +  Dac  +  Da^c  +  k'c  _kc 

dr2  r  3r  3z2  1  A  2 


d3*C  +  D3C  +  d^C 
3r2  r  3r  3z2 


+  k  -  k  C  -  0 
1  2 


[40] 


Equation  [40]  is  subject  to  the  boundary  condition 


C  -  C*  -  kiA2.  r  -  z  -  <*>,  all  t. 


[41] 


and  we  again  consider  initially  the  simplest  condition  at  the  electrode 
surface 


C  -  0,  0  <  r  <  a,  z  -  0,  all  t 


[42] 


and  with 


D 


ac 


dz 


0,  r  >  a,  z  -  0,  all  t 


[43] 
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The  similarity  of  Equation  [40]  in  standard  form 


d2C  f  1  dC  t  a2c 
BtZ  x  dx  dz2 


[44] 


is  apparent  and  we  obtain  the  solution  for  the  concentration  of  B  as 


C(r,z,k2) 


JQ(ar)ada 


f 


Q(p,k2)JQ(ap)pdp 


[45] 


By  analogy  to  the  distribution  of  sinks  for  the  chronoamperometric  transient 
we  assume  that  for  the  c.e.  reaction  we  can  write 


-Q(«) - ± - 

a  cos (0  ) 


d  (k  )  +  d  (k  )cos(9)  + 
0  2  12 


+  d  (k  )cos"(0) 

n  2 


[46] 


I 


where  the  d  are  now  functions  of  k  .  Then  at  z  -  0  we  obtain 
J  2 


A 
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n 


J-0 


Var)Jfj< 

\ (aa)  d(aa) 

n 

H 

[47 

[k  1 1/2 

2  2  ,  . i ( J-l > 12) 

a  +  5  J  <aa) 


Then,  if  the  boundary  condition  [42]  applies,  the  (n+1)  coefficients  d  (k  ) 

J  2 

can  be  derived  by  solving  the  set  of  equations 


r  148 ' 

2 


r"  i  J  (ar  )  J  /  \  i 

J ‘  (?) 


(aa)  d(aa) 


<  < J-X ) /2  ) 


at  (n+1)  radial  positions.  A  working  curve  of  the  kinetically  limited  flux 
can  then  be  constructed  from 


V  vv 

R(k  )  -  2*a  )  — 

Z_-  <J+1> 

j-0 


[49] 


Mass  transfer  of  the  reactant  species  F  in  the  e.c.  (catalytic)  reaction 


scheme 
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E  +  ze 


F  +  G 


(iii) 


E  +  H 


Is  described  by 


D  i!c  +  2  ac  +  D  .  k'c  c 

dr2  r  dr  dz2  2  G 


ar2  r  ar  dz2 


kC-0 

2 


in  the  presence  of  excess  of  the  substrate  G.  Provided  the  reaction  (iii)  can 
be  driven  sufficiently  hard  in  the  forward  direction,  [50]  will  be  subject  to 
the  boundary  conditions 


C-C  ,  0<r<a,  z-0,  all  t 


C  -  0,  r  -  ®,  z  -  ®,  all  t 


where  C°°  is  the  concentration  of  species  E  in  the  bulk  of  the  solution.  With 
a  distribution  of  sources  over  the  disk  having  the  form  [46] ,  we  therefore 
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obtain  the  concentration 


C(r,z,k2) 


■if 


exp 


r  2  k2i 

1,2. 

L 

r® 


P*  ?] 


k  T  1/2 


Jq  (ar)ada 


Q(p.k2)JQ(ap)pdp 


[53; 


and  at  z  -  0 


r-  f  Jo(ar)Jfjli^a 

-  -  )  VV  _ L.d. . . 

Da^  Jo  r  2  \y'l  , 

j-o  a  +  n  (<»a) 


(aa)  d(aa) 


[54] 


((J-D/2) 


The  d^(k2>  are  evaluated  by  solving  this  set  of  equations  and  the  working 
curve  of  the  kinetically  limited  flux  can  again  be  constructed  using  149], 

The  use  of  the  boundary  conditions  [42]  and  [51]  is,  however,  unrealistic 
for  the  same  reasons  that  [14]  is  unrealistic  for  the  evaluation  of  the 
chronoamperometric  transients.  In  reality,  the  electron  transfer  reactions 
(ii)  and  (iii)  will  take  place  at  finite  rates  and  the  e.c.  (catalytic) 
reactions  especially  will  be  subject  to  the  effects  of  the  distribution  of 
potential  in  the  solution.  It  is  better  therefore  to  seek  to  predict  the 
complete  form  of  the  polarization  plots  taking  into  account  all  these  effects. 
It  is,  in  fact,  possible  to  achieve  such  a  complete  description  of  the  systems 


21 


since  we  are  considering  their  behavior  in  the  steady  state;  we  restrict 
attention  here  to  the  case  where  the  electron  transfer  processes  are 
irreversible  so  that  they  can  be  described  by  the  boundary  condition  [10]. 

The  distribution  of  potential  in  the  solution  follows  as  outlined  above 
and  we  obtain  at  z  -  0 


n 


J-O 


[55] 


Application  of  the  boundary  conditions  [10]  to  the  reaction  step  (ii)  in  the 
c.e.  mechanism  taking  into  account  [55]  gives 


'  dz  J  a  cos(0) 


d  (k  )  +  d  (k  )cos(0)  + 
0  2  12 


1 


+  d  (k  )cosn(0) 

n  2 


-  k  exp 


fa*<r,k2)F] 


RT 
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L 


-  k  exp 


-zaF  a 
xRT 


d  (k  )e  f  - 
J  2  i  J la 


J-o 


i  i  J  (ar) J t  \ 

Izl 

Da£-  Ja  r2  Vl1'2 


(aa)  d(oa) 


j-0 


r  2  2]  1  .  .  (  (J-l)/2)  I 

L°  +  D  J  (aa> 


[56] 


where  k  is  che  rate  constant  for  reaction  (iii)  for  zero  ohmic  potential  at 

the  surface  of  the  disk.  The  non-linear  equation  system  [56]  can  be  solved  at 

(n+1)  radial  positions  for  the  (n+1)  required  coefficients  d  in  the  same  way 

J 

as  has  been  described  for  the  case  of  the  tertiary  current  distribution  in 

(1  2) 

simple  redox  reactions  at  disk  electrodes  '  . 

In  a  similar  way  the  application  of  [10]  and  [55]  to  the  reaction  step 
(iii)  of  the  e.c.  (catalytic)  mechanism  gives  the  equation  system 


a  cos (8) 


d  (k  )  +  d  (k  )cos (.6)  +  •••  +  d  (k  )cos  ( 8 ) 
0212  n  2 


-  k  exp 


\  d  (k)ef  II 

*RT  ZlJ  2  J  J^a 

J-o 
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j-o 


J0(«r)J 


M 


( aa )  d(aa) 


r  k  t  1/2  I 

2.  2  .  . ( ( j- 1 ) / 2)  I 

+  D  J  (aa)  ' 


[57] 


where  k  is  now  the  rate  constant  for  reaction  (iii)  (assumed  here  to  be 
irreversible)  for  zero  ohmic  overpotential  in  the  solution.  The  coefficients 
derived  from  [56]  and  [57]  are  naturally  also  functions  of  k  (i.e.  of  the 
standard  rate  constant,  of  a  and  of  the  overpotential)  as  well  as  of  the 
solution  conductivity  k.  The  polarization  curves  can  therefore  be  completely 
described  for  this  model  (and  indeed  other  models)  using 


R(k2,k,a) 


dj  (kj ,k»*) 

(J+D 


[58] 


Blsgtfs.slgp 

It  can  be  seen  that  the  application  of  Neumann's  integral  of  two 
parameters  (the  radial  and  s-dependence  of  the  surface  source/sink  over  the 
disk  electrode)  allows  a  more  comprehensive  discussion  of  the 
chronoamperometric  transients  than  has  been  achieved  to  date  using  other 
analytical  as  well  as  simulation  techniques1*’7’9  16 5  .  In  particular,  since 
Neumann's  integral  leads  to  a  relation  between  the  assumed  form  of  the  flux 
and  the  concentration  oat  the  surface  (Equations  [17]  and  [18])  it  becomes 
possible  to  derive  the  transient  for  boundary  conditions  such  as  [10]  which 
are  more  realistic  than  the  assumption  of  zero  or  constant  surface 
concentration,  [14],  which  has  been  used  to  date.;  for  example,  Equation  [2 
leads  to  the  definition  of  the  potentiostatic  relaxation  behavior  of  disk 
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electrodes . 

We  have  pointed  out  elsewhere'*'8,16'  that  the  combined  effects  of  the 
finite  rates  of  electrode  reactions  and  of  the  potential  distribution  in  the 
solution  must  make  the  distribution  of  the  flux  over  the  surface  much  more 
uniform  than  that  predicted  using  the  boundary  condition  [14].  Indeed,  the 
predictions  of  a  model  of  uniform  surface  flux  and  zero  average  surface 
concentration<4,16>has  been  shown  to  be  In  close  accord  with  experimental 
measurements<7> .  While  the  methodology  outlined  here  allows  a  systematic 
exploration  of  one  of  the  effects,  that  of  the  finite  rates  of  electrode 
reactions,  Equation  [10],  the  exploration  of  the  effects  of  changes  In  the 
potential  in  the  solution  unfortunately  remains  restricted  to  the  discussion 
of  the  relaxation  behavior<17>(cf.  Equations  [34]  and  [35]);  in  effect,  we 
reach  the  limits  of  analytical  techniques.  However,  it  does  become  possible 
to  develop  a  complete  discussion  of  linear  sweep  voltammetry  and  of  cyclic 
voltammetry(17) ,  as  has  been  outlined  for  the  first  case  above,  provided  the 
effects  of  the  tertiary  current  distribution  are  neglected. 

It  has  been  shown  previously  that  there  are  no  restrictions  in  the 

(1  2 ) 

discussion  of  the  steady  state  behavior  of  simple  electrode  processes 
The  discussion  presented  here  shows  that  this  is  equally  true  for  electrode 
processes  coupled  to  reactions  in  solution  such  as  for  the  c.e.  and  e.c 
(catalytic)  mechanisms;  a  complete  description  of  such  systems  can  be 
developed.  Comparisons  of  the  predictions  with  those  based  on  other  models 
and  methods  of  analysis<*’18,19'2°>as  well  as  extensions  to  other  reaction 
schemes  will  be  presented  elsewhere'215 . 
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glossary  si  Symhala 


a 

b 

1 

c 

1 

d 

i 

e 

i 

f 

l 

C 

C°° 

c 

Av 

c’ 

D 

F 

i 

o 

i 

I 

J 

l 

k 


k 

k 


Q 

R 

r 


T 


Disk  radius,  cm 

Weighting  function  series  coefficients 

Fourier  coefficients,  mols  (cm  s)  1 

Fourier  coefficients,  mols  (cm  s)  1 

Fourier  constant  terms 

Fourier  series  integral  terms,  cm*1 

Concentration,  mols  cm  3 

Bulk  concentration,  mols  cm*3 

Average  concentration,  mols  cm*3 

Surface  concentration,  mols  cm*3 

2  *1 

Diffusion  coefficient,  cm  s 

Faraday  constant,  96485  C  equivalent’1 

.  j 

Exchange  current  density,  A  cm 
Current,  A 

*2 

Current  density,  A  cm 
Bessel  functions 

Heterogeneous  rate  constant,  cm  s  1 
Homogeneous  rate  constant,  various 
Heterogeneous  standard  rate  constant,  cm  s 
Homogeneous  rate  constants,  various  types 

Flux,  mols  cm  2s  1 
Gas  constant,  8.314  J  mols  lK  1 
Radial  coordinate,  cm 
Laplace  transform  variable 
Temperature,  K 
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z  Coordinate  normal  to  plane  of  disk,  cm 

z  Charge  of  an  ion 

o  Transfer  coefficient  (when  in  exponent) 

a  Continuous  dummy  integration  variable . 

Solution  potential,  V 
p  a  sin(0) ,  cm 

k  Solution  Conductivity,  ohms  ^m  1 

ij  Overpotential ,  V 

r  Laplace  transform  dummy  integration  variable 
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